Tunneling in Fractional Quantum Hall line junctions 
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We study the tunneling current between two counterpropagating edge modes described by chiral 
Luttinger liquids when the tunneling takes place along an extended region. We compute this current 
perturbatively by using a tunnel Hamiltonian. Our results apply to the case of a pair of different 
two-dimensional electron gases in the fractional quantum Hall regime separated by a barrier, e. g. 
electron tunneling. We also discuss the case of strong interactions between the edges, leading to 
nonuniversal exponents even in the case of integer quantum Hall edges. In addition to the expected 
nonlinearities due to the Luttinger properties of the edges, there are additional interference patterns 
due to the finite length of the barrier. 

PACS numbers: 73.43.-f, 71.35. Lk, 71.23.An 



INTRODUCTION 



Edge states in the fractional quantum Hall effect 
(FQHE) are very interesting examples of one-dimensional 
strongly interacting quantum systems. The right and left 
moving edge modes of a quantum Hall bar are spatially 
separated by a macroscopic length and this leads to ex- 
ponentially small backscattering. The main reason for 
localization is thus suppressed and the edge modes are 
a nearly ideal ballistic system. Many experiments have 
been devoted to the study of their unique characteris- 
tics pj . An interesting geometry is a constriction of the 
electron gas, the so-called quantum point contact. By 
means of an electrostatic potential created by a gate, 
the two edges of a bar are brought in close proximity, 
allowing tunneling phenomena to take place at a single 
point of the fluid. Recently,progress in the technique of 
cleaved-edge overgrowth 0, Lj has led to the fabrication 
of samples in which the tunneling now occurs along a bar- 
rier of mesoscopic extent between two spatially separated 
two-dimensional electron gases. Kang et al. 0, El have 
performed detailed studies of the conductance of these 
new structures. Their samples consist of two-dimensional 
electron gases (2DEGs) separated by an atomically pre- 
cise barrier of length 100 /im and of width 8.8 nm. They 
have studied the conductance of this structure as a func- 
tion of the applied magnetic field and the voltage bias be- 
tween the two gases. Many theoretical works have tried 
to explain their results Hi,0lll![I3[HEEi0- 

In this paper, we give the results of a perturbative 
calculation of tunneling between two edge modes that 
are counterpropagating. Each of these edge modes are 
described by a chiral Luttinger liquid and we focus on 
the situation where they are characterized by the same 
anomalous exponent g. We consider the situation where 
tunneling takes place along an extended region with con- 
stant amplitude. The anomalous exponent g of the chiral 
Luttinger liquid is governed by the bulk FQHE fluid(s) 
and enter in the expression of the correlation function of 
the particle that tunnel. The extended tunneling geome- 
try is potentially rich of new interference phenomena not 
found in single point-contact devices. Notably, T. L. Ho 



pointed out the existence of oscillating currents without 
AC drive in the case of integer edge modes [l^. We find 
that in the FQHE regime there is a nontrivial interplay 
between the well-known nonlinearities of the Luttinger 
liquid and the interference pattern of the tunnel barrier. 
We also discuss the case of strong interactions between 
the left and right-moving modes (but still weak tunnel- 
ing) which is amenable to the same theoretical treatment. 
The exponent g then takes nonuniversal values dictated 
by the details of the interaction potential, even for in- 
teger quantum Hall edges. The formulas we obtain are 
generalizations of previous results already available in the 
literature |16U17| . This geometry may also give rise in the 
presence of disorder to a derealization transition [TsL Il9| . 
The line junction geometry may also be relevant to a re- 
cent set of experiments using samples with an extended 
constriction [2(|. A detailed modelling of such line junc- 
tions has been performed by Papa and McDonald [2lll22| . 

In section II, we introduce a model of edge states in the 
quantum Hall regime with an extended tunnel barrier. In 
section III, we present the perturbative treatment of the 
weak tunneling regime. Section IV gives the results of 
our study. Finally section V contains our conclusions. 



II. THE GEOMETRY AND KINEMATICS OF 
TUNNELING 

The geometry we consider is illustrated in fig. ^ It is 
that of an extended barrier whose height remains con- 
stant along some spatial extent we call L. The barrier 
separates spatially two 2DEGs. This is the geometry of 
the samples studied by Kang et al. 0]. In the quantum 
Hall regime there are edge modes of conduction that are 
counterpropagating. With the atomically precise barrier 
of ref. |2|, there is tunneling only in the integer quan- 
tum Hall regime when the applied bias between the two 
2DEGs is small. This geometry may also be realized by 
electrostatic gates, in which case the energetics may be 
different and extended tunneling may eventually be real- 
ized also in the fractional quantum Hall regime. 

We first discuss the kinematics of tunneling by reason- 
ing in the case of bulk filling factor v = 1 . This situation 
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was first investigated by T. L. Ho [15|- The Landau lev- 
els are degenerate in the bulk and their energies raise 
when they approach a barrier. This scheme applies to 
both sides of the barrier. Ultimately the Landau levels 
of both sides will cross in the forbidden region inside the 
barrier and tunnel effect will lead to the opening of single- 
particle gaps. A simple effective theory clearly displays 
these phenomena. We first introduce the left and right 
moving chiral fermions if>R and ipL that are the relevant 
modes when the bulk filling factor is between v = 1 and 
v = 2. They are the counter-propagating edge states [l6| . 
Their kinetic energy purely due to the confining potential 
is then : 



dx(ili R d x tp R - ip{d x ip L ), (1) 



where v is the drift velocity along the barrier. In mo- 
mentum space the left modes have energy ei(fc) = —vk 
and the right modes have e R {k) = vk. These dispersion 
relations cross at k = for zero energy. It is at this point 
that the tunnel effect is strongest. The tunneling through 
the extended barrier can be described by a Hamiltonian 
mixing these modes : 



^-tunnel 



+L/2 



L/2 



dx (tp R (x)ip L (x) +ipl(x)ip R (x)), 



where T is the tunnel amplitude, constant along the bar- 
rier by hypothesis. This is the tunneling Hamiltonian we 
treat in this paper. The coordinate x is defined as in 
fig- HI so that the modes propagate for all values of x and 
tunneling is restricted to the range (— L/2, +L/2). 
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Figure 1: Two laterally separated 2DEGs under a magnetic 
field. A tunnel barrier causes mixing of the counterpropa- 
gating edge modes. Realistic configuration (a), topologically 
equivalent geometry (b) 

This Hamiltonian can be seen a mass term of the 
fcrmion field for an infinite-length barrier. This mass 
leads to a gap in the free fcrmion excitation spectrum 
and this gap implies the existence of a range of energies 
for which there are no propagating solutions but only 
evanescent waves along the barrier. This phenomenon 



repeats in a real sample also at the crossings of the other 
higher- lying Landau levels as discussed by Kang et al.. 

The problem defined by 'H; ,,, ■ 'H, „,,„,! is a free fermion 
theory and has been studied exactly [l5( • The scattering 
properties of the barrier can be deduced by a straightfor- 
ward Landauer-typc calculation. However this simplicity 
does not extend into the realm of the FQHE : here the 
electron operators are described by a chiral Luttinger liq- 
uid (xLL) theory and no longer by the free theory Eq.JTJ. 
We study the case with no interactions through the bar- 
rier so the xLLs are decoupled and all the properties of 
their correlation functions are known. We still use the 
tunneling Hamiltonian Eq.J2J and treat it in perturba- 
tion to obtain the current - voltage I — V characteristic 
of this system. We consider the case of a single edge ap- 
propriate to a v = 1/m, m odd, FQHE fluid. A single 
chiral boson is then enough to describe the X-LL theory. 
The bulk fluid determines the value of the anomalous 
exponent g characterizing the correlations at the edge. 
When the tunneling is between two distinct gases then 
the operators ip R and ipL should be taken as electron 
operators. Contrary to the case of the quantum point 
contact geometry, the bosonized version of the theory 
cannot be treated exactly since the tunneling term is no 
longer a boundary operator. We are thus left only with 
the perturbative calculation. 

In real samples, there may be also noticeable Coulomb 
interactions along and across the barrier which can be 
modeled by the following kind of Hamiltonian : 



Hint = J dxdy Vi{x - y)[p R (x)p R (y) + p L (x)p L {y)] 
+ 2 / dxdy V 2 (x - y) p R (x)p L (y), (3) 



where p R (x) =: ip R (x)ip R (x) : is the right-moving edge 
electron density operator describing density fluctuation 
at point x (pl(x) is similarly defined). This may be po- 
tentially important in the structure of Kang et al where 
the width of the barrier is smaller than the magnetic 
length. In this case, Wen ^(| has shown that the eigen- 
modes are a mixture of the uncoupled left and right mov- 
ing edge modes. He showed also that the study of the 
tunnel effect may be performed along the standard lines 
with simply a redefinition of the Luttinger parameter g 
when the momentum dependence of the interactions may 
be neglected. So the results of our perturbative tunneling 
calculation also apply to interacting edges if we use an a_ 
propriate value of g. For example, Mitra and Girvin 
have estimated the value ofo to be 0.6 — 0.7 in the 
structures studied in refs.(@, 111)- This value is an ex- 
ample of what can be expected in the integer quantum 
Hall regime when interactions destroy the perfect quan- 
tization of g by the bulk physics. Strictly speaking, this 
analysis applies when the interedge is translation invari- 
ant along the barrier, like in Eq.@ and the tunnel 

We note that such perturbative calculations have been 
done in various tunneling situations and are a valuable 
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tool to perform spectroscopy of reduced dimensionality 
samples]!! IH El |H E3 in the quantum Hall regime. 

In a given sample under a magnetic field, the edge 
states will be populated up to a definite Fermi level. We 
will consider the symmetric situation where the Fermi 
level is the same on both sides of the sample. Applica- 
tion of a bias voltage induces a difference between these 
Fermi levels. The external magnetic field may be tuned 
so that the Fermi points meet at k = 0. This is the case 
first observed by Kang et al where there is a zero-bias 
peak in the conductance. The peculiar dispersion of the 
edge modes leads to special kinematic constraints in the 
case of an extended barrier. Indeed, when the length of 
the barrier goes to infinity, momentum along the barrier 
is a conserved quantity during tunneling (contrary to the 
case of point tunneling). Also energy conservation in the 
tunnel effect means that only states at the intersection 
of the dispersion relations of the left and right moving 
modes will tunnel. For a nonzero Fermi wavevector (de- 
fined with respect to the "vacuum" situation Eas.i|ll2|)) 
this implies that one needs a finite bias eV to get a tun- 
neling current. The voltage threshold will go to zero if 
the filling factor is tuned to the degeneracy point kj — 0. 



III. PERTURB ATIVE TREATMENT 

We now describe the application of the standard tun- 
neling formalism to the line junction. Tunneling occurs 
from a non-equilibrium situation in which the chemical 
potential of the left edge fiL is different from fi R of the 
right side because of the applied voltage eV — [II — MR- 
The current is expressed as the rate of change of the num- 
ber of particles Nl — J dx(ip L (x)ipL(x)), for example, on 
the left-hand side of the sample. Thus, the average tun- 
neling current is given by: 



I(t) = -e(N L (t)), 



(4) 



where e is the electron charge. We set Ti = 1 and ks = 1 
everywhere. 

The total Hamiltonian is of the form : 



ti = H.l + tij 



ti 



tunnel j 



(5) 



where H.l and 7i R are Hamiltonians for the chiral Lut- 
tinger liquids on each side of the junction. We treat 
titunnei as a perturbation and keep only the leading term. 
Since we limit ourself to the tunneling at low voltage, low 
temperature and small tunneling amplitude T, titunnei 
is the only term which does not commute with Nl as 
Hl and Tin separately conserve Nl and Nr. Standard 
first order expansion in the interaction representation 
yields H| : ' 



Following |28|, this procedure leads to the following 
formula for the tunneling current : 



2 , dkdk' sm 2 [{k~k')L/2] ^ 



4tt z 



(k - k') 



l\2 



(7) 



de 
2^ 



A R {k,e)A L (k',e - eV){n f (e - eV) - n f (e)}. 



In this equation n/(e) = 1/(1 + exp(/3e)) is the Fermi 
factor at temperature T {(3 = 1/T) and A R (k,oj) (resp. 
Ai J (k,uj)) is the chiral spectral function for the right 
(resp. left) moving chiral Luttingcr liquid. 

The spectral function may be obtained from the imag- 
inary part of the retarded Green's function : 



AR,L(k,U}) 



2 Im 



dtdx 
4tt 2 



Atjjt— ikx 



(8) 



-i0(*)Wfl,L(M),< L (O.O)}>/3 



where we take the thermal average (}p in a grand canoni- 
cal ensemble including chemical potentials (J-l.r- We first 
evaluate the Fourier transform in space and time of the 
forward Green function L of right-movers and left- 
movers; then the spectral densities are obtained by use 
of the following identity which holds for fermions at finite 
temperature : 



G> L (k, u) = (1 - n f (u))A RtL (k, u). 



(9) 



The Green functions for a chiral Luttinger liquid are 
given by : 



G>{k,u>)= a^ 1 



dtdx 
4tt 2 



D iu>t—ikx±ikfX 



[i sinh(7rT(i ^x/v))) 9 '' 



(10) 



where g is the Luttinger liquid parameter ,16j, kf is the 
fermi wavevector and a is a microscopic cut-off. In the 
case of electron tunneling, the anomalous exponent is 
g = 1/v = m when the two FQHE fluids have filling 
v = 1/m. This will describe the situation where a barrier 
separates two distinct electron gases. Here we consider 
only the simplest situation with the same filling factor on 
both sides of the barrier. If we consider tunneling with 
interactions through the barrier, then g is non-quantized, 
e. g. in an integer Hall system it may be slightly smaller 
than one (due to repulsive interactions). 

The Fourier transforms of the spectral functions can 
be computed (2^| : 

G> L (k, lo) = § (uj + v(k f T k)) 



I(t)=ie / dt'{[N L {t),Htunnel(t')]). 



(6) 



(11) 
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This leads to the following closed formula : 
A R , L (k, u) = 2(— y- 1 cosh(-/3w) 5{u + v(k fT k)) 



XiJi — h z , % , 

K 2 2ir ' 2 2tt j ' 



(12) 



where B(x,y) is the Euler Beta function. In the zero 
temperature limit, formula 1|12|) reduces to the standard 
zero-temperature chiral Luttinger liquid spectral func- 
tion : 



(k, uj) = a 



9-1 



2vr 



kTkf^^Siuj + vikfTk)). (13) 
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Figure 2: The scaling function A g (z)of a chiral Luttinger liq- 
uid as a function of the scaling variable z = u/T in the case 
8=0.7. 

The spectral function can reexpressed to show its scal- 
ing properties (we define the Fermi energy e / = vkf) : 




T 



A R<L (k,u) = 2(2tt) 9 - 1 Ij-j A g (u/T) 8{u+e f ^vk), 

with the scaling function : 
A g (z) = (afc,)*" 1 cosh(z/2) B{ 9 - + i ± | - i ±). (15) 

We note that the u/T scaling is destroyed when j ^ 1 
by the overall factor (T / e f) g ~ l . The scaling function A g 
is plotted in Fig. [5] 



where we have defined an overall current scale Iq = 
eT 2 (ak f ) 29 - 2 /{%TT£ f T(g) 2 ). It is not possible to compare 
directly currents with and without interactions because 
the current scale Jo is a function of g. More precisely 
the microscopic cut-off appears through the dimension- 
less combination akf when g ^ 1. In the g = 1 case, one 
can perform the remaining integral in Eq. l|lt)|) in terms 
of the Sine integral [30( : 



io 



k f L [Si(2k f L{v + 1)) + Si{2k f L(v - 1))] (17) 
1 - cos(2fc/L(¥+l)) 1 - cos(2fc /j L(tJ- 1)) 



2(v + l) 



2(17-1) 



A. Non-resonant case 

We first discuss the case when kf ^ 0. This means 
that, in the absence of bias, either the two Fermi points 
are both below the maximum tunneling point k = 
where the dispersion relations do cross or both are above 
the k — 0. So tunneling is suppressed till a finite bias. 
This can be seen in formula Eci. (|16fl . If we set the barrier 
length to infinity then the sine function in the integrand 
of the tunnel current Eq. I|16(l peaks to a delta function. 
The current is then nonzero only beyond a bias threshold 
given by \eV\ > 2vkf. This is simply due to the conser- 
vation of momentum in the limit L — > oo which severely 
restricts tunneling. This threshold is smoothed out for 
finite L values and disappears in the point contact limit 
L — > 0. We will refer to the generic case kf ^ as being 
"non-resonant". On the contrary, for kf — tunneling is 
allowed for infinitesimal bias, the "resonant" case treated 
in the next section. In a realistic set-up one needs a spe- 
cial fine-tuning of the magnetic field to reach the reso- 
nant case. In the experiment of Kang et al, there is an 
extended region of kf, i.e. of magnetic field where there 
is sizeable tunneling for infinitesimal bias as revealed by 
a zero-bias peak in the conductance. Such an extended 
range of tunneling is out of reach of the present model 
where there is a single point satisfying the resonance con- 
dition. 

In fact, the tunnel current is O(L ) below the threshold 
and grows as O(L) above. Below the threshold, the L — > 
oo limit leads to : 



IV. RESULTS 

Our formula Eq.© yields the value of the tunneling 
current for any barrier length L, interaction parameter 
g, finite applied voltage and temperature. It is convenient 
to define a dimensionless voltage v = eV/(2sf). We will 
concentrate on the zero temperature limit. The tunnel 
current can be written as : 



/ = Jo (k f L) 



du 



sin 2 [(l + u)kfL] 
(\+u) 2 (k f L) 2 



2\S-1 



(16) 



io 



2g-l 

2e s ) 2r( 5 + l/2) 



eV\' y ~ l y/*T{g) 



,3 1 feV 
x 2 F l( l,-, ff+ -^ — 



(18) 

where 2J1 is the hypergeometric function |3(l |. Above 
threshold, we find the simple formula : 



(19) 



Some J — V curves are displayed in fig. (|3J) for a very 
long barrier, kfL^l. In the FQHE case, we have taken 
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Figure 3: Tunneling current vs eV/sf for three values of g : 
g = 3, (dotted line) g — 1 (solid line), <? = 0.7 (dashed line) 
for a barrier of length kfL — 1000. Since the normalization 
Io of the currents is a function of g, no comparison can made 
between the absolute value of these I — V curves. 



the Luttingcr parameter g = 3 for tunneling between two 
v = 1/3 liquids. The case g = 1 refers to free fermions, 
here there is a saturation of the tunnel current, i. e. I 
becomes O(L) independent of V. Finally, we have drawn 
the I —V curve for a Luttinger parameter less than one, 
possibly relevant to the case of IQHE edge modes with 
interedge interactions. There is a sharp current spike 
at threshold. The current indeed behaves as ~ (eV — 
2ef) 9 ~ 1 close to threshold in the large- L limit. For V — > 
0, the current behaves as V 29 ^ 1 . 

For a finite length barrier with kfL of order unity, the 
kinematic singularities are smoothed out. Typical curves 
are presented in Fig.@. The behavior for V small is 
unchanged but there are additional oscillations due to 
the diffraction pattern generated by the barrier acting 
like a slit. We note that in the case with g less than one, 
the singular behavior for small voltage oc V 29 " 1 shows 
up only for small voltage. In fig.JSJ is presented a close- 
up of the case g = 0.7 at very small bias : here one 
can observe the nonlinearity typical of a Luttinger liquid 
and the additional interference pattern. The period of 
the oscillations is 2ir/(k/L) in terms of the dimensionless 
voltage eV/(2e f ). 

Fig.Elshows the tunneling current vs the length of the 
barrier for a fixed voltage eV/(2sf) — 0.1 and several 
values of the interaction parameter g. There is a satura- 
tion for L larger than hundreds of the Fermi wavelength 
A/, as discussed above, because the voltage bias is below 
the threshold. There is also a beating pattern in both 
cases. The fast oscillations are due to the interference of 
L versus the Fermi wavelength while the slow oscillations 
are ruled by the voltage scale introduced in the problem 
by the bias voltage eV/(2ef) set to 0.1. We have omitted 
the case g = 3 because it has a similar shape but with a 
very small current in units of Iq. 




eV/£ f 

Figure 4: Tunneling current vs eV/ef for three values of g : 
g — 0.7, g — 1 and g = 3 from top to bottom for a barrier of 
length k f L= 10. 
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Figure 5: Tunneling current vs eV/ef for the interacting pa- 
rameter g = 0.7 for a finite barrier length kfL = 1000. 
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Figure 6: Tunneling current in unit of 7o vs length of the 
barrier for the interacting parameter g = 1 (solid line) and 
g = 0.7 (dotted line). The bias is fixed at eV/(2e f ) = 0.1. 



Figure 7: The auxiliary function governing the crossover be- 
tween the regimes of long and short barriers for values of the 
parameter g. 



B. Resonant case 

We now turn to the case kf = which can be ob- 
tained by tuning the external applied magnetic field. 
This means that in the absence of bias the left and right 
Fermi levels exactly coincide at k = 0. The tunnel cur- 
rent from our perturbative calculation is now simpler be- 
cause the Fermi energy disappears from the problem. As 
a consequence we find : 



eT 2 L 2 



aeV 



2g-l 



f g (- 



,eVL, 



T(g) 2 8irav \ 2v J JHV 2v 
where we define the auxiliary function f g (x) by : 



+ i 



dt (i - r) 



\xt) 



(xtf 



(20) 



(21) 



This function decreases as ir/x for large x and any value 
of g and has a nonzero value at the origin. This means 
that I cx V 29 " 1 if eVL/v < 1 while I cx V 29 " 2 if 
eVL/v 3> 1. For a typical g = 0.7 value this means that 
the current should first rise as V 0A before going down 
at larger voltage as V~ 6 . The function that governs 
the crossover is shown in Fig.JJJ). Some I — V curves are 
displayed in Fig.©. While the Luttinger liquid nonlin- 
earities are still present, the oscillations formerly due to 
the presence of the scale kf no longer exist. The g = 0.7 
curve shows that strong interedge interactions in a situ- 
ation of weak tunneling lead to a diverging conductance 
at zero bias. 

Finally we note for completeness that our study may 
be extended at nonzero temperature by using the finite- 
temperature spectral functions. Using a dimensionlcss 
temperature z = we have : 



eT 2 (ak f ) 
8ire f 



du ■ 



r(.g) 2 



-i 2 - 

z 

)k f L] 



2g-2 



r(.# 

4tt 2 



(kfL) 



(1 



'(kfL)' 



e zu (l + e' ziu+v ^){l + 



(22) 

,-z{u-v)\ 
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Figure 8: Tunneling current vs bias eVL/(2v) at resonance 
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We expect thermal rounding of the oscillatory features as 
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well as the threshold behavior in the non-resonant case. 



V. CONCLUSIONS 

In this paper we have studied the tunneling between 
two counterpropagating edges modes pertaining to differ- 
ent electron gases by means of a perturbative expansion. 
We find the characteristic nonlinearities of the Luttinger 
liquid properties. There are additional oscillations due to 
the finite extent of the barrier. Our study applies to the 
case of tunneling between two FQHE fluids for v = 1/ m 
and also applies to the case of tunneling between edges 
when there is strong interactions throughout the barrier 
(but weak tunneling) in which case the characteristic ex- 



ponent entering the spectral functions is no longer quan- 
tized by the bulk physics. It remains to be seen if samples 
with the correct properties can be artificially produced. 
We note that the present experiment by Kang et al is 
apparently in the regime of strong tunneling (even if the 
conductance throughout the barrier is far below the ex- 
pected Landauer value). 
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